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We explore charge transfer in the telomere G-Quadruplex (TG4) DNA theoretically by 
the nonequilibrium Green’s function method, and reveal the topological effect of charge 
transport in TG4 DNA. The consecutive TG4(CTG4) is semiconducting with 0.2 
0.3ey energy gap. Charges transfers favorably in the consecutive TG4, but are trapped 
in the non-consecutive TG4 (NCTG4). The global conductance is inversely proportional 
to the local conductance for NCTG4. The topological structure transition from NCTG4 
to CTG4 induces abruptly ~ 3nA charge current, which provide a microscopic clue 
to understand the telomerase activated or inhibited by TG4. Our findings reveal the 
fundamental property of charge transfer in TG4 and its relationship with the topological 
structure of TG4. 
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1. Introduction 

Charge transfer along double-stranded DNA has attracted much attm tion among 
biomedicine, chemistry and physics communities in the past decadeP^ As a special¬ 
ized DNA sequence, telomeres involve many essential physiol ogic al processes, DNA 
damage, cell replication, aging, genetic stability and cancers^ Human telomeres 
consist of tandem repeats of the hexanucleotide {TTAGGG)n 5-lOkb in length 5’-3’ 
toward the chromosome end, terminating in a single-stranded 3’-overhang of 100- 
200 basest^ In normal cells, the telomere shortens with each cell replication. When 
a critical length is reached, cells undergo apoptosis.l^ Nevertheless, 80-90% of can¬ 
cer cells preserve their telomere length by the activation of telomerase and thus 
become abnormal.^ Interestingly, the G-rich sequence can fold into G-quadruplex 
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(G4) (see in Fig.l), which is a secondary structure consisting of stacked G-tetrad 
planes connected by a network of Hoogsteen hydrogen bonds and stabilized by centre 
monovalent cations, such as Na~^ and K~^. For human telomeric DNA, the repeats 
of the telomere (TTAGGG)n form the telomere G4 (TG4), where intramolecular 
G4 forms three-layer packet (Fig. 1),^ which can form several topological structures. 
It has been found that the formation of TG4 inhibits telomerase activity, which ob¬ 
structs tumor immortal mechanism.l^ Hence, TG4 has attracted extensive studies 
as an attractive target for cancer therapeutic intervention.l^l 

Moreover, because of DNA self-assembly properties, its paring specificity and 
conformational flexibility offer great potential for the rational design of DNA-based 
nanostructure and nano electronics Different topological structures of G4 DNA 
provide more possibilities for nano devices. 

However, there are divergent opinons on charge transfer in the G-rich DNA se¬ 
quences. Sugiyama et. al. believe that the G base is favor for charge transfer because 
the ionization potential of G base is lowest among the nucleobases.miThe numerical 
study also shows G4 DNA favoring charge transport based on a simplified mono-G 
tight-binding G4 model.n^ In fact, charge transfer in G4 could be quite different 
from normal DNA due to the cation effect as a new hydrogen bond in G4 DNA. 
Barton’s group found experimentally that G4 has great trapping potency on charges 
by chemical florescence method.!!^ Actually, the problem of charge transfer in G4 is 
still open because of the complicated topological structures of G4. The full under¬ 
standing of charge transfer in G4 requires comparing different topological structures 
of G4 DNA. However, it is not easy to compare functions of different topological 
structures of G4 by biomedical methods. Physical method can provide a microscopic 
insight to understand the mechanism of the G4 functions in physiological process 
and charge transfer for nano devices. 

In this paper we will study charge transfer in intramolecular TG4 by the effective 
tight-binding model with nonequilibrium Green’s function (NEGF) method. Inves¬ 
tigating the current-voltage (I-V) characteristics of two predominant topological 
structures of TG4 in vivo, we find novel charge transfer properties of TG4 beyond 
the classical charge transfer law, which provide basic properties for designing G4- 
based nano devices and inspire a physical mechanism on the telomerase activation 
inhibited by TGiW 


2. General G4 model 

The TG4 DNA can be viewed as comprising of four parallel tt stacks for conducting 
charge channels through the superposition of tt orbitals along DNA molecules, which 
is observed by the NMR spectroscopy and x-ray crystallography.tSI We consider the 
two ends of G4 DNA connected with two semi-infinite one-dimensional electrodes. 
The Hamiltonian of G4-DNA model can be written as 


I 


I 


H — HOi -b Hl,R -b HCi-LR + Henv 


( 1 ) 
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The charge channels in the backbone of G4 DNA is very small and the tt orbit 
channel dominates charge transfer in G4 DNaIISI Thus, the effective tight-binding 
Hamiltonian of G4 DNA can be expressed as 


N 


^ ^ ^n,n+l yC^Cn-\-l T 

n—l,k£G4 

— y] tm (fkCl + cjcfj 
(k,l)&Gi 


( 2 ) 


where Snk = £« + 4 ^^ on-site energy of each basespH where the second 

term is the ion effect in G4 packet, where i?o = 1.99nm is the distance between 
the ion and tt electrons, and Sr ~ 2eo is the effective dielectric constantP^ The 
second term in Eq.dS]) describes the tt electron hoppings between Gs in G4_packet. 
The tn^n+i is the nearest tt electron hopping parameters listed in table iJHUHl xhe 
notation 5’-XY-3’ indicates the direction along the DNA strand (see,Fig.1(b) 

The tm describes the hopping amplitude between Gs in G4 packet. The Hlh in 
Eq.([T]) is the electrodes at the left and right ends of G4 DNA, 


Hlm 


'y ^ tc Cn-l-1 + 

n^L,R 


( 3 ) 


The Hg^-lr in Eq. o is the couplings between the G4 DNA and electrodes, 

Hg^-lr = -th (cjci -I- cjco) - tR ^cj^CAT+i -I- . (4) 

The c^(cfc) is the creation (annihilation) operator of electron or hole at k sites. 
The £l(r) ( tL{R) ) is the on-site energy (hopping parameter) in the electrode, 
respectively. To minimize the contact effect, we use the optimal injecti on c ondition 
\/tatL{R) = the strong coupling in the electrode tL(_R) = leE,^^ where a 

denotes A, C, G, and T. 


Table 1. The on-site energy and the tt electron hop- 


ping parameter between nearest bases. 

5’-XY-3’(eV) 

X/Y 

G 

C 

A 

T 

G 

0.119 

-0.295 

-0.186 

0.334 

C 

0.026 

0.042 

-0.008 

-0.161 

A 

-0.013 

0.091 

-0.038 

-0.157 

T 

0.044 

-0.066 

-0.068 

0.180 

Sn 

7.75 

8.87 

8.24 

9.14 


Since realistic DNA molecules are under physiological conditions, we consider 
TGi DNA having water and ion environment. These water and ion environments of 
TG 4 provides a stochastic electric field background to influence the TG4 structure 
as a stochastic perturbation. We introduce the Hgnv describing these environment 
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effects by the stochastic fluctuation of the G 4 structure. Thus, the Hgnv has the 
same form to Hci, but the parameters are e®"’' = 0.5Ar and = 0.5tn,n+i^r, 

where A describes the stochastic strength and r € (—1,1) randomlyl ^^ l ^^ l 

It should be remarked that we introduce the electrode such that we can calculate 
the TV characteristic, which can be compared directly with experiments. For the 
case of G4 DNA in vivo, suppose that the physiological environment provides an 
electric potential bias to G4 DNA. Thus, the situation of G4 DNA in vivo maps 
similarly to the case we consider. 

Using NEGF Method, the current can be expressed in terms of^^H 

^ ^ - /«) (5) 

where G’’(“^ = [E — Hg^ — the retard (advanced) Green’s 

Function; the = *(Xl(_r) “Xl(_r)) describes the level-width function, where 

XL(fl) is the retard (advanced) self-energy of electrodes. The Jl^r) is the Fermi 
function of the left (right) electrodes. 

3. TG4s and their classifications 

In general, the telomere sequence can be folded to many types of TG4 DNA. The 
human telomere {TTAGGG)n has been found to be folded to two topological struc¬ 
tures (TPl and TP2) in vivo shown in Fig.l.^SlThe TG4 packet as an basic element 
can be constructed to different configurations of TG4 DNA, which can be classified 
into three classes, the consecutive mono(hybrid)-TG4(GM(H)TG4) DNA, and non- 
consecutive mono(hybrid) TG4 (NGM(H)TG4) DNA. Moreover, we may compare 
the TG4 packet with and without ion. For convenience, we label the TG4 DNA by 
sn{ijk)i or sn(ijk)a, where sn labels the telomere sequence (TTAGGG)n\ the (ijk) 
labels the TG4 topological structure and the last i and a labels ion and absence 
of ion. For example we define sl 6 = TTA{GGGTTA)iq and when sl 6 is folded to 
two TPls and TP2s with ion. We label it by sl6(1122)i that is GHTG4. When sI 6 
is folded to one TPl and two TP2 without ion. We label it as sl6(1022)a that is 
NGHTG4, where 0 labels the sequence of non-TG4 packet. 

4. Topological effect of charge transfer 

In order to capture the basic physics of the charge transfer in TG4, we first turn 
off the environment effects Henv = 0. For GMTG4, s4 = TTA{GGGTTA)^, it can 
be folded to only one TG4 packet. We compare the TV characteristic of the s4 
telomere chain and two topological structures (TPl and TP2) of TG4 with and 
without ion, namely s4(l)i, s4(2)j, s4(l)a, and s4(2)a, shown in Fig.2 (a). The s4 
chain is of little conductance compared with CMTG4. The s4(l)i(a) and s4(2)i(a) 
show semiconducting with the energy gap 0.3eV. The currents saturate at around IV 
due to only finite states for molecule systems. The saturated currents of s4(l)i and 
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Fig. 1. (Color online) The sketch of two topological structures of TG4 DNA in human vivo. 


s4(2)i reach 6nA and InA respectively, which is much higher than that of s4(l)a, 
and s4(2)o (0.06nA). That means that the ion in TG4 enhances charge transfer. 
For longer telomere sequences, sl2 = TTA{GGGTTA) 12 , and its GMTG4, the basic 
behaviors of the I-V characteristics are similar to that of the s4 case, but the energy 
gaps for the non-ion cases are very small. The saturated currents of s4(2)a is larger 
than that of s4(l)a, but reverse for the sl2 case, which can been seen in Fig.2 (b). 




V (V) V (V) 


Fig. 2. (Color online) The I-V characteristic of CMTG4 DNA in (a) for the s4(j)f(a), in (b) for 
the sl2{jjj)i{a), where j = 1,2 means TPl or TP2. 


For CHTG4 the basic physical behavior of I-V characteristics is same to that of 
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CMTG4. They are also semiconductor with ^ 0.3eV energy gap and the saturated 
currents of the case with ion are much larger than that of the cases without ion. 
The saturated currents depend on the configurations of the topological structures 
of TG4 shown in Fig.3. 
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Fig. 3. (Color online) The I-V characteristic of CHTG4 DNA in (a) for the sl2(jkj)i(a), and in 
(b) for the where j = 1,2 means TPl or TP2. 


For NCHTG4 in Fig.4, the saturated currents of NGHTG4 with and without 
ion become very small 0.002 ^ 0.2pT and actually vanish for some cases, which is 
quite different from GMTG4 and CHTG4. The saturated current of the telomere 
chain is larger than that of NGM(H)TG4. It implies that TG4 in NCM(H)TG4 can 
trap charges, which is consistent with the experimental results.ll^ 
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Fig. 4. (Color online) The I-V characteristic of NCMTG4 DNA in (a) for the sl2{jj0)i{a) and 
sl2(0j0)i(a), and in (b) for the sl6{jjk0)i{a), sl6{0jj0)i{a) and sl 6 ( 00 j 0 ) 2 (a). 


Gharge transport in NCM(H)TG4 exhibits novel topological properties. In com- 
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mon sense of charge transport in classical and quantum systems, the global con¬ 
ductance of the system is proportional to the local conductance of the system. 
However, for NCH(H)TG4, the global conductance is inversely proportional to the 
local conductance. The conductance of TG4-DNA system may be regarded as to be 
proportional to the saturated current. We find the conductance Gtga > GsA from 
results in Figs.2 and 3. When the telomere chain sl6 forms NGMTG4, sl6(1000)f, 
the conductance becomes Gsi6(iooo)i < shown in Fig.4. It implies that the 
global conductance of NGM(H)TG4 is inversely proportional to the local conduc¬ 
tance even though the local conductance Gtga > GsA- In other words, charge 
transfer in NCM(H)TG4 exhibits nonlinearity. The topology of TG4 induces the 
anomalous charge transport. We can understand these novel charge transfer prop¬ 
erties in NCM(H)TG4 from quantum mechanics. Physically, mobile charges form 
Bloch wave in a periodic or quasi-periodic structures of CM(H)TG4, while mobile 
charges in NCM(H)TG4 are trapped due to Anderson’s localization effect in disor¬ 
der chain. That is why charges transfer favorably in the CM(H)TG4. Actually, DNA 
sequence may be regarded as a quantum system. Gharge transfer in DNA sequence 
has been found exhibiting some quantum nature.l^ The anomalous conductance in 
NGM(H)TG4 we find is a novel quantum property of charge transfer induced by 
topology of TG4 DNA. 

5. Topological structure transition and mechanism on telomerase 
activation and inhibition 

The charge transfer properties of TG4 DNA can give a physical clue to understand 
TG4 how to activate or inhibit telomerase. It has been found that 80% ^ 90% can¬ 
cer cells obtain their immortality by activating telomerase to preserve the telomere 
length.l^The telomerase activation or inhibition is related t o TG4 DNA.^^The chal¬ 
lenging problem is how to activate or inhibit the telomerase !^^!^^! Suppose that the 
ligands of the telomerase offer an electric bias on the telomere DNA to form a circuit 
inducing charge current. The topological structure transition from NGM(H)TG4 to 
GM(H)TG4 generates abruptly the charge current 3nA from 0.002pA (see Figs. 
2 and 4). This current as a threshold current can activate or inhibit telomease. 
This scenario can provide a clue to understand TG4 blockading the cancer cell im¬ 
mortal pathways and eventually perishing tumors in the biomedical experimental 
observation.!^ This finding inspires a further challenging problem for biomedicine 
how to realize or control the topological structure transition between NGM(H)TG4 
and CM(H)TG4. 

6. Discussion 

It should be remarked that the ratio of the saturated currents between GM(H)TG4 
and NGM(H)TG4 has 6 orders. The solvent leads only O.I ~ 0.5eV on-site energy 
fluctuation to DNA bases, which cannot lead to higher G4 structures.!^ We in¬ 
vestigate the environment effects in Fig. 5 for two cases. When the environment 
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Fig. 5. (Color online) The ratio of the currents of CM and NCM TG4s at V = IV. We take 20 
samples to give the average and the standard errors. The values of the ratios still remain 6-order 
differences against the environment effects. 


fluctuation is small A < 0.5, the error bar of the ratio of the currents for the 
S'12(lll)i/S'12(110)i is larger than that of the S12(222)i/S12(202)i, which implies 
that the current fluctuation depends on the detail configurations of the TG4 in the 
small environment fluctuation. When the environment fluctuation becomes larger, 
A > 0.6, the current fluctuation trends to be insensitive to the detail configurations 
of TG4. More importantly, we do see that the 6-order current difference induced 
by the topological structure transition from NGM to CM TG4s is robust against 
the environment fluctuation. Moreover, it has been found that two basic mecha¬ 
nisms of charge t ran sfer in DNA chain, coherent tunneling (hopping), and thermal 
diffusive hopping .1231 The coherent tunneling or hopping occurs betwe n G /C pairs 
and the thermal hopping occurs in more than 3 consecutive A/T pairs.l23lFor TG4, 
{TTAGGG)n, the coherent hopping dominates the charge transfer due to only 3 
consecutive A/T pairs in TG4. Thus, the coherent transport method is valid for 
TG4. In addition, the energy gap at Fermi energy is about 0.2 ~ 0.3eV for TG4, 
which is I order higher than the influence of temperature fluctuation in room tem¬ 
perature, and the charge transfer in TG4 is not sensitive to temperature fluctuation 
because the coherent hopping dominates the charge transfer in TG4 instead of the 
thermal diffusive hopping. Namely, temperature fluctuation in room temperature is 
not sensitive for our results. Therefore, in room temperature all our conclusions are 
invariant. 
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7. Conclusion 

In summary, we reveal novel charge transport properties in TG4 DNA. The 
GM(H)TG4 is semiconducting with 0.2 ^ 0.3eV energy gap. The conductances of 
different configurations of TG4 DNA follow the inequalites, GcM(H)TGi > Gpc > 
GNCM(H)TGi) where TG means telomere chain. The conductances depend on the 
topological structures and configurations of TG4. The cation in TG4 enhances 
charge transfer. The TG4 packet in NCM(H)TG4 suppresses charge transfer, which 
agrees with the experimental results.^I^ The topological structure transition be¬ 
tween NCM(H)TG4 and CM(H)TG4 induces a threshold charge current to activate 
or inhibit telomease. Our findings not only reveal novel charge transfer properties 
in TG4 DNA that can offer many opportunities for DNA-based electronics, but 
also provide physical insights and hints to understand the telomerase activation 
inhibited by TG4. 

The biomedical investigation found that tumor growth is related to the telom¬ 
erase activation and inhibition. The telomerase activation and inhibition is related 
to TG4.It implies strongly that states of TG4 could be related tumor growth 
even cancer. However, the biomedical investigation cannot find the working mech¬ 
anism of these three factors, TG4, tumor growth and telomerase activation and 
inhibition. Our results provide an understanding of the working mechanism of the 
telomerase activation and inhibition. We found that the topological structure tran¬ 
sition from NGM(H)TG4 to CM(H)TG4 generates abruptly the charge current 3nA 
from 0.002pA (6 order difference) when the ligands of the telomerase offer an elec¬ 
tric bias on the telomere DNA. It strongly implies that this abrupt current as a 
threshold current can activate or inhibit telomease. Namely, if we can control the 
topological structure transition from NGM(H)TG4 to GM(H)TG4 we can control 
the tumor growth. This scenario can provide a hint to understand TG4 blockading 
the cancer cell immortal pathways and eventually perishing tumors in the biomed¬ 
ical experimental observation.^^ Gonsequently, the further challenging issue is how 
to tune the topological structure transition of TG4 by biomedical method, which is 
expected to further study. 

These novel charge transfer properties in TG4 DNA give some fundamental re¬ 
lationships between topology, anomalous charge transport and biomedical function. 
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